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Abstract 

We study some dynamical properties of a Lorentz gas. We have considered both the static and 
time dependent boundary. For the static case we have shown that the system has a chaotic com- 
ponent characterized with a positive Lyapunov Exponent. For the time-dependent perturbation 
we describe the model using a four-dimensional nonlinear map. The behaviour of the average 
velocity is considered in two situations (i) non-dissipative and (ii) dissipative. Our results show 
that the unlimited energy growth is observed for the non-dissipative case. However, when dis- 
sipation, via damping coefficients, is introduced the senary changes and the unlimited engergy 
growth is suppressed. The behaviour of the average velocity is described using scaling approach. 
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1. Introduction 

The process in which a classical particle acquies unlimited energy growth from collisions 
with heavy boundaries, also called as phenomenon of Fermi acceleration, were first reported by 
Enrico Fermi QJJ] as an attempt to explain the acceleration of cosmic rays. He proposed that such 
behaviour was due to interaction between charged particles and time-dependent magnetic fields. 
Soon after some alternative models have been proposed using different approaches with ap- 
plication in different fields of science including molecular physics [2], optics [3], nanostructures 
33], quantum dots fl. 
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One of the most studied versions of the problem is the well known one-dimensional Fermi- 
Ulam model (FUM) J2, 0, B H [HI Qjl Q2] • The model consists of a classical particle confined 
and bouncing between two rigid walls in which one of them is fixed and the other one moves 
according to a periodic function. It is well known that the phase space, in the absence of dis- 
sipation, shows a mixed structure in the sense that depending on the combinations of control 
parameters and initial conditions, both invariant spanning curves (also called as invariant tori), 
chaotic seas and Kolmogorov- Arnold-Moser (KAM) islands are all observed. An alternative 
model was latter proposed by Pustylnikov 11131 1 1 4fl - Such system consists of a classical parti- 
cle bouncing in a vertical moving platform under the effect of an external constant gravitational 
field 115, la, 17, 18, lH 2(1 21]. Both models seems to be quite similar. However there are 
many differences between them. The main difference is, that in the FUM framework the Fermi 
acceleration is not observed. On the other hand, for specific combinations of both control param- 
eters and initial conditions the phenomenon of unlimited energy growth can be observed in the 
bouncer model. This apparent contradictory result was latter discussed and explained by Licht- 
enberg and Lieberman |22, 23tl an d can be easily understood by looking at the space phase. The 
FUM has a set of invariant spanning curves limiting the size of the chaotic sea (as well as the 
particle's velocity), but such invariant tori, which could be interpreted as a barrier, they are not 
observed in the bouncer model and the energy grows unbounded. However, in two-dimensional 
systems it is not so simple to say if the phenomenon of Fermi acceleration will be observed or 
not. In this sense a conjecture was proposed by Loskutov-Ryabov-Akinshin (LRA) 12411 . Such 
conjecture, known as LRA-conjecture, states that a chaotic componet in the phase space with 
static boundary is a sufficient condition to observe Fermi acceleration when a perturbation is in- 
troduced. Results that corroborate the validity of this conjecture include the time dependent oval 
billiard JH], stadium billiard JH]. 

When dissipation is introduce, it causes a drastic change on the phase space. Invariant span- 
ning curves are destroyed, the elliptic fixed points turns into sinks and the chaotic sea is replaced 
by a chaotic attractor [27]. However, the influence of dissipation on the averave velocity is still 
not fully understood. It was discussed by Leonel lp"oll some results for a one-dimensional Fermi- 
Ulam model. It is well known that such system in its original formulation the particle do not have 
unlimited energy growth. If one consider the motion of the moving wall to be random the phe- 
nomenon of Fermi acceleration is observed. However, the introduction of inelastic collision is 
enough to suppress Fermi acceleration. Results considering also the one dimension Fermi-Ulam 
model under an external force of type sawtooth IB ill under effects of dissipation generated from 
a sliding of a body against a rough surface has also been considered. The external perturbation 
of sawtooth type was chosen because the oscillating wall always gives energy to the particle after 
collisions. The main question is: Would it be possibel to suppress Fermi acceleration under the 
efect of dissipation generated from a sliding of a body against a rough surface? The answer is not 
so simple and it dependes on both, initial conditionand the combination of control parameters. It 
was observed that for a certain range of parameters Fermi acceleration indeed happen, however, 
it is suppresses under specific conditions. 

In this paper, we study the problem of a Lorentz Gas considering both the static and the 
time-dependent boundary. In the first part of our work we study the Lorentz Gas with static 
boundary. We derived a two dimensional nonlinear mapping that describes the dynamics of 
the model. We obtain the phase space and we show that it is chaotic with positive Lyapunov 
Exponent. Then, in the second part of our paper we introduce a time-dependent perturbation on 
the boundary. There are many different ways to introduce time dependent perturbation and the 
most common methods are: stochastic case, where the boundary changes according to a random 
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function l28ll ; regular case, where the size of the boundary varies according to an harmonic law 
lf29ll . However, in both situations the center of mass it is assumed to be fixed. Now, for the 
first time, we introduce a different kind of time dependent perturbation for a Lorentz gas. We 
assume that the radius of the scatters are fixed and the center of mass changes according to an 
harmonic function. Our main goal in this part of our work is to verify the validity of the LRA 
conjecture, which is confirmed when we studied the behaviour of the average velocity for an 
ensemble of particles. Since the phenomenon of Fermi Acceleration is present in this model our 
next approach is to introduce dissipation into the model via damping coefficients and trying to 
understand what is the influence of dissipation on the particles behaviour. Our results allow us to 
confirm that when inelastic collisions is introduced into the model it is a sufficient condition to 
break down the phenomenon of Fermi acceleration. In both, conservative as well as dissipative 
case, we describes the behaviour of average velocity using scaling formalism. 

The paper is organized as follows. In section|2]we describe how to obtain the two-dimensional 
mapping that describes the dynamics of the static system. Section[3]is devoted to discuss the time 
dependent model as well as our numerical results. Finally, conclusion and acknowledgments are 
drawn in section [4] 



2. A static Lorentz Gas and the mapping 

In this section we discuss all the details needed for the construction of a non-linear mapping 
that describes the dynamics of the problem. The model consists of a classical particle of mass m 
suffering elastic collisions with circular scatters (see Fig[T]i. We choose a triangular arrangement 
of the scatters (which also seems to be the Star of David if one connect circles 1,5 and 9 and 1 1,7 
and 3) in order to avoid particles traveling infinitely far between collision and the fixed lattice 
spacing a to be twice the radius of the scatters, R-a/2. The system is described in terms of a 
two dimensional mapping E(8„,b„) = (6 n+ \,b n+ \) where the dynamical variable 6„ denotes the 
direction of the trajectory while b n is the impact parameter. Given an initial condition (0Q,bo), 
the particle starts from the black circle (center in Fig. [TJ and hit one of the 12 others circles. In 
this sense, we specify the scattered hit in the collision n+1 by s„ = 0, 11 and we introduce 
l{s„) for the distance of this scatterer and the scatterer hit at the collision n+1. l(s„) can assume 
the values of 2a/ V3 and 2a for even and odd values of s„, respectively. Additionally, when the 
particle hits the boundary it is specularly reflected with the same absolute velocity. The particle 
does not suffer influences of any external field along its linear trajectory. From the green triangle 
in Fig|2](a) one can easily verify 



. ,„ ns n (b„ + i -b„) 

sin(<9„ - — ) = — — . (1) 

6 l(s„) 



Moreover, from Fig. |2](b) one can find that 



a = aicsin(-b n+ i/R), (2) 

P = n-2a, (3) 

iff = 2a -6 n , (4) 

9 n +i - n-2a + 6„. (5) 



Such a result allow us to obtain 9„ +i , 




O n+ i - n + 6„ + 2 arcsin(^„+i//?). (6) 
From Eq. [TJit is easily to find that the impact parameter, b n +\, is given by 

ns n 

b„+\ =b„- l(s n )sm(0„ - —). (7) 
o 

Thus, the mapping that describe the dynamics of a two dimensional Lorentz gas is given by 

0„+i = n + 9„ + 2 arcsin(%^) 

b n+1 = ft B - l(s n )sm(6 n - ^) ' 1 ; 

by definition & e [-i?, /?] and 9 e [0, 27r] is a counterclockwise angle such that S is defined on 
the fundamental domain x [0,2^]. From the mapping S, Eq. ([8), one can easily obtain 

the Jacobian Matrix, J, which is defined as 



ae„ +i de„+i \ 
ae„ db„ 

db„ +l 0b„ +l 
de„ db„ 



(9) 



with coefficients given by 

d6 n+l 
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= l-2^ 2 -^ +1 Z(s„)cos(0„--^), (10) 
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^ 2 -b 2 n+l , (ID 
^ = -/(,„) C os(0„-^), (12) 
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Figure 2: Dependence of (a) b n+ \ on b„ and 8„; (b) 8„ + i on b, 1+ i and B n . 



After some easy calculation one can show that the mapping E preserves the phase space 
measure since det(7) = 1 . 

It is well know the Lyapunov exponents are an important tool to identify whether the model 



is chaotic or not. As discussed in [32], the Lyapunov exponents are defined as 

Aj = lim - In | A ,-| , ; = 1,2 , (14) 

«-»oo n 

where Aj are the eigenvalues of M = Y[1=i Ji(fih bi) an d Ji is the Jacobian matrix evaluated over 
the orbit (0,-, £>,). However, a direct implementation of a computational algorithm to evaluate Eq. 
(fT4l has a severe limitation to obtain M. Even in the limit of short n, the components of M 
can assume very different orders of magnitude for chaotic orbits and periodic attractors yielding 
impracticable the implementation of the algorithm. In order to avoid such problem we note that J 
can be written as J = &T where is an orthogonal matrix and T is a right triangular matrix. Thus 
we rewrite Mas M — J„J n -\ ■ ■ .Ji®\®i l J\, where T\ = ®~ l J\. A product of J2&1 defines anew 
J' T In a next step, it is easy to show that M = J n J n -\ ■ ■ ■ J^®!®^ 1 J^Ti- The same procedure can 
be used to obtain T2 = & 2 l J' 2 and so on. Using this procedure the problem is reduced to evaluate 
the diagonal elements of Tj : T' n , T' 2r Finally, the Lyapunov exponents are now given by 

1 " 

Aj= Urn -XlnirU , j= 1,2 . (15) 

i=\ 

If at least one of the Aj is positive then the orbit is classified as chaotic. Additionally, in conserva- 
tive systems A \ + A2 = and in dissipative systems Ai + A2 < 0. The phase space for the mapping 
© is shown in Fig. [3 a). For the same control parameters used in Fig. [3](a), a = 2, we have also 
evaluated numerically the positive Lyapunov exponent as one can see in Fig. [3b). The average 
of the positive Lyapunov exponent for the ensemble of the 5 time series gives A = 0.371 ± 0.001 
where the value 0.001 corresponds to the standard deviation of the five samples. 



3. A Time Dependent Lorentz Gas. 

In this section we introduce a new kind of time dependent perturbation in a Lorentz gas. We 
assume that the radius of the scatters are fixed and the center of mass changes according to an 
harmonic function /(f), in particular we assume the case where 



/;(*) = e,[l + cos(f)], for i = x,y. 
5 



(16) 
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Figure 3: (a) Phase space generated from iteration of mapping {8); (b) Behaviour of the positive Lyapunov exponent of 
the chaotic sea. The control parameter used in both figures were a = 2. 



where e,- is the amplitude of the time-dependent perturbation and t is the time. When we intro- 
duce time dependent perturbation into the model, two new dynamical variables appears, namely, 
velocity and time. Now, we have a two dimensional system described in terms of a four dimen- 
sion nonlinear mapping which relates the collision n with the (n + Y) ,h , i.e., H(6>„, b n , V„, t„) = 
{6„ + i,b„ + \,V„+\,t„ + i). The corresponding variables are: the direction of the trajectory, n ; the 
impact parameter, b„; the absolute velocity of the particle, V n and the instant of the hit with the 
boundary, t„. 

Assuming that an initial condition (f?o, bo, Vo, ?o) is given, we can obtain the equation that de- 
scribe the dynamics of the system. Thus, according to our construction, the cartesian components 
of R are given by 

X(6 n ,t n ) = Rco&(S n ) + e x [l +cos(*„)L (17) 
Y(6„,t n ) = Rsm(6 n ) + e y [l +cos(f„)]. (18) 

where S„ is the angular position which is given by 6„ = n/2 + 8„- aicsm(b„/R). Since we already 
know the angle that the particle's trajectory does with the horizontal (9„ + nil) and the position 
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of the hit at the collision n , we can obtain the vector velocity of the particle that is written as 

f „ = |V^|[cos(0„ + 7r/2)T+ sin(0„ + n/2)J\ , (19) 

where i and j represent the unit vectors with respect to the X and Y axis, respectively. The above 
expressions allow us to obtain the position of the particle as a function of time for t > t n : 

X p (t) = X(6 n ,t n ) + \f n \cos(0 n +n/2)(t-t n ), (20) 
Y p (t) = Y(6 n ,t n ) + \P n \sm(e n +n/2)(t-t n ). (21) 

The index p denotes the corresponding coordinates of the particle. In order to know the 
position of the particle at (n + \) th collision we need to solve numerically the following equation 



r = ^j[X x (t) - X p (t)] 2 + [7,(0 - Y p (t)f = R , (22) 

where both X x and Y y are given by 

X x (t) = l x + e x [l + cos(0], (23) 
Y y (t) = l y + e y [l + cos(0] , (24) 

being l x and l y the X and Y components of l(s„); this distance is measured from the origin of the 
coordinates system to the center of the s n = 0, 11 scatters at (n + \)' h collision. Since we al- 
ready know the position of the particle at the collision (n+l) th , one can easily find the distance be- 
tween two successive impacts, which is given by d = y/[X p (t) - X(6„, f„)] 2 + [Y p (t) - Y(6 n , f„)] 2 . 
Then, the time at (n + \)th collision is obtained evaluating the expression 



VLMO - X(6 n , t n )Y + [Y p (t) - Y(S n , t n )Y 

tn+l =t„-\ ^ . (25) 

|V„| 

The next step is to obtain the impact parameter, b n+ \, which is given by 

b n+ \ =b n - l(s n )sm {e n -iff+ | j, (26) 

where l(s„) = y/(AX) 2 + (AY) 2 and if/ = arctan(AX/AT) with 

AX = l x + e x [cos(t n+1 ) - cos(f„)] (27) 
AY = l Y + e,,[cos(f„ + i) - cos(f„)] (28) 



Moreover, the new direction of the trajectory, 6 n+ \ is 
6 n +\ = n + 6 n + 2 arcsin 



R 



(29) 



We already know (6 n +\,b n+ \,t n+ \), however we still have to find"^„ + i. At the new angular 
position 5„+i, the unitary tangent and normal vectors are 

7*„ +1 = cos(5„ + i); + sm(S n+ i)j , (30) 
7^ B+ i = -sin(£ n+ i)z' + cos(£ n+ i),/'. (31) 



Since the referential frame of the boundary is moving, then, at the instant of the collision, 
according to our construction, the following conditions must be matched 

V*„ + i-7*„+i = yV' n -^ n+1 , (32) 
V' n+1 -lt„ +1 = -8V' n -Tt n+1 , (33) 

where ye [0, 1] and S e [0, 1] are damping coefficients, which means that the particle can lose 
velocity/energy upon collision. The complete inelastic case occurs when y = 6 = 0. On the other 
hand, when y = S = 1 corresponds to the conservative case. The upper prime indicates that the 
velocity of the particle is measured with respect to the moving boundary referential frame. 
Hence, one can easily find that 

= y^n-^n + l+(l-y)^b(tn + l)-^n + l ■ (34) 



t„ +1 • Tt n+1 = -SP„ ■ ~N„ + i + (1 + 8Y$b{t n+ i) ■ ~$ n+ i , (35) 
where Vt(t n +i) is the velocity of the boundary which is written as 

~$b(t„+i) = - sin(f„ +1 )[e c / + e y j] , (36) 
Finally, the velocity at (n + \)' h collision is given by 

«+i •7* n+ i) 2 + (t„ +I -^„ +1 ) 2 . (37) 

3.1. Numerical Results 

Our numerical results for the time-dependent Lorentz Gas shows basically the behaviour of 
the average velocity of the particle. Two different procedures were applied in order to obtain 
the average velocity. Firstly, we evaluate the average velocity over the orbit for a single initial 
condition which is defined as 

y -^TlZ y - (38) 

where the index i corresponds to a sample of an ensemble of initial conditions. Hence, the 
average velocity is written as 

1 M 

V =M^ V - (39) 

where M denotes the number of different initial conditions. We have considered M = 1000 in 
our simulations and from now on, we also fixed the value a = 2. 
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Figure 4: Behaviour ofVxn for different initial velocities. The control parameters used were e x = 1(T 4 , €y = 3 X 10 
and a = 2 



3.2. Scaling results for the conservative case 

Our main goal in this section is describe a scaling present in the model for the conservative 
case, where, in Eq. [32] and Eq. [33] we assume y — 5 — 1, and also verify the validity of LRA 
conjecture. 

We begin discussing a scaling observed for the average velocity of the particle as function of 
Vo and n. It is shown in Fig. [4] the behavior of the V x n for different initial velocities. Hence, 
the control parameters used in Fig. [4]are e x = 10~ 4 , e y = 3 X 10~ 4 . Additionally, 21 different 
values of Vo were chosen and for each one we randomly chose t e [0, 2tt], 9 e [Q,2n] and 
b G [1 - (e x + €y), -1 + (e t + e^)]. As one can see, all curves of the V behave quite similarly 
in the sense that: for short n, the average velocity remains constant, then after a changeover, all 
the curves start growing with the same exponent. Such behaviour is typical in systems that can 
be described using scaling approach. Based on the behavior shown in Fig. [4] we propose the 
following hypotheses: 

1 . When n <k n x , V behaves according to 

V sar oc V ( , (40) 
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2. For n » n x , the average velocity is given by 

V ocn v , 



(41) 



3. The crossover iteration number that marks the change from constant velocity to growth is 
written as 

n x (x v* , (42) 

where £ and v are the critical exponents and £ is a dynamic exponent. 

After consider these three initial suppositions, we suppose that the average velocity is de- 
scribed in terms of a generalized homogeneous function of the type 

V(V ,n) = lV(FVo,l' ? «), (43) 

where 1 is the scaling factor, p and q are scaling exponents that in principle must be related to (, 
v and If we chose Wo = 1, then 1 = V l ^ p and Eq. d43l is given by 

V(.V ,n) = V - 1 ' p V 1 (.\r q/p n), (44) 

where V\(V^ q ^ p n) = V{\, VZ q ^ p n) is assumed to be constant for n <k n x . Comparing Eq. ( l44b and 
Eq. ( I40l i. we obtain £ = -I /p. 

On the other hand, Choosing now 1 = n Eq. d43l is rewritten as 

V(V , n) = n- 1/9 y 2 (n _p/9 Vo) , (45) 

where the function V2 is defined as V2(n~ p ^Vo) = V(n~ p ^ 9 Vo, 1). It is also assumed to be constant 
for n » n x . Comparing Eq. (05]) and Eq. (|4TT > we find v = -1/g. Given the two different 
expressions of the scaling factor 1, we obtain a relation for the dynamic exponent which is 
given by 

(46) 



Note that the scaling exponents are determined if the critical exponents £ and v were numer- 
ically obtained. The exponent v is obtained from a power law fitting for the average velocity 
when n » n x . Thus, an average of these values gives v = 0.49(1). Figure [5] shows the behaviour 
of (a), V sat X Vo and (b), n x x Vo. Applying power law fittings we obtain f = 1.00(1) = 1 and 
£, = 2.01(3). Considering the previous values of both ( and v and using £ = f/v, we find that 
£ = 2.04(2). Such result indeed agrees with our numerical data. In order to confirm the initial 
hypotheses and, since the values of the scaling exponents v and £ are now known, we will 
collapse all the curves onto a single and universal plot, as demonstrated in Fig. [6] Additionally, 
such result allow us to confirm the validity of LRA conjecture since the V grows unbounded for 
n » n x . 

3.3. Scaling results for the dissipative case. 

In this section we will use the same scaling formalism as used in the previous section. We 
concentrate to characterize the behaviour of the average velocity in terms of the number of colli- 
sions with the boundary and as a function of the damping coefficient along the normal component 
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of the particle's velocity, 6. We study a dissipative version of the Lorentz Gas close to the transi- 
tion from unlimited to limited energy growth. Indeed, such a transition happens when the control 
parameter 6 — > 1 and it is better characterized if adopt the following transformation 8 — > (1 — 6). 
To obtain the average velocity, each initial condition has a fixed initial velocity, Vq = 10~ 4 and 
randomly chose t e [0, 2tt], e [0, 2n] and b e [1 -(e x + e y ), -1 +(e x + e y )]. The control parameter 
y were fixed as been y — 1 . 

It is shown in Fig. Qthe behaviour of average velocity as function of the number of collision 
for different values of 6. Note that, for different values of 6, the average velocity, for small 
«, starts to grow with the same slope and them they bend towards a regime of saturation for 
long enough values of n. The changeover from growth to the saturation is marked by a typical 
crossover number n x . For such a behaviour, we can propose the following scaling hypotheses: 

1 . When n <k n x the average velocity is 

V oc n\ (47) 



2. For long time, n » n x , the average velocity approaches a regime of saturation, that is 
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Figure 6: (a) Behaviour of average velocity for different values ofVo; (b) their collapse onto a single and universal plot. 

described as 

V M oc(l-Sf, (48) 

3. The crossover number that marks the regime of growth to the constant velocity is written 

as 

n x ^(\-8f, (49) 

where cr, rj and z are critical exponents. 

These scaling hypotheses allow us to describe the average velocity in terms of a scaling 
function of the type 



V[n, (1 - 6)] = lV[l p n, 1 9 (1 - 6)] , (50) 

where p e q are scaling exponents and 1 is a scaling factor. Since 1 is a scaling factor, we can 
chose it such that V'n - 1 , yielding 

V[n,(l-S)] = n- l,p Vi[(ri)- q/p (l-S)] , (51) 
12 
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Figure 7: Behaviour ofVx nfor different values of 5, as labeled in the figure. 



where Vi[(ri)~ q ' p (l — d)] = V[l,(ri)~ q ' p (l—5)] is assumed to be constant for n <K n x . Comparing 
Eqs. (l47l i and (f5TJ, we obtain 77 = -1/p. A power law fitting gives us that 77 = 0.471(2). Such 
value was obtained from the range of S e [0.99, 0.99999]. Choosing now - 8) = 1, we have 
that 1 = (1 — SY l l q and Eq. (l50l l is rewritten as 

V[n, (1 - $)] = (1 - <5r 1/9 y 2 [(l - , (52) 

where V2[(l - S)~ p ' q n] = V[{1 - 5)~ p ' q n, 1] is assumed to be constant for n » n x . Comparing 
Eqs. (1481 ) and d52l . we obtain -l/q — cr — -0.47(3) [see Fig. |8](a)]. Using now the expressions 
obtained for the scaling factor I, we can easily show that 

z = - = 0.99(5) , (53) 

which is quite close to the value obtained numerically, as shown in Fig. [8](b). A confirmation of 
the initial hypotheses is made by a collapse all the curves of V x n onto a single and universal plot, 
as shown in Fig. [9] showing that the system is scaling invariant. We also shown that dissipation 
causes a drastic change in the behavior of V. Note that when S — > 1, implies that Eq. [48] and 
Eq. [49] diverge, thus recovering the results for the conservative case, i.e., Fermi acceleration. 
However, when 6 is slightly less than 1, the average velocity grows and then reach a regime of 
saturation for long enough time. Our results reinforce that dissipation introduced via damping 
coefficients is a sufficient condition to suppress the phenomenon of Fermi acceleration. 



4. Conclusion 

In this paper we consider the problem of a classical Lorentz Gas considering both the static 
and time-dependent boundary. For the static case we obtain the mapping that describes the dy- 
namics of the system and we have shown that the model has a chaotic component characterized 
with positive Lyapunov Exponent. After that, we introduce a new type of the time-dependent 
perturbation on the boundary. Our results confirm the validity of LRA conjecture for the con- 
servative case since the phenomenon of Fermi acceleration is observed. When dissipation, via 
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damping coefficients, is introduced into the model we observed that the average velocity grows 
with time and then reaches a constant value for large enough time confirming that Fermi accel- 
eration is suppressed. Finally, the average quantities were described by scaling functions with 
characteristic exponents whose validity was confirmed with the collapse of the curves of V into 
an universal plot. 
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